Quantum Opt. 2 (1990) 253-265. 



Generation of discrete superpositions of coherent 
states in the anharmonic oscillator model 

A Miranowicz, R Tanas and S Kielich 

Nonlinear Optics Division, Institute of Physics, Adam Mickiewicz University, 
Grunwaldzka 6, 60-780 Poznan, Poland 

Abstract. The problem of generating discrete superpositions of coherent states in 
the process of light propagation through a nonlinear Kerr medium, which is modelled 
by the anharmonic oscillator, is discussed. It is shown that under an appropriate choice 
of the length (time) of the medium the superpositions with both even and odd numbers 
of coherent states can appear. Analytical formulae for such superpositions with a few 
components are given explicitly. General rules governing the process of generating 
discrete superpositions of coherent states are also given. The maximum number of 
well distinguished states that can be obtained for a given number of initial photons is 
estimated. The quasiprobability distribution Q(a, a* , t) representing the superposition 
states is illustrated graphically, showing regular structures when the component states 
are well separated. 



Quantum Opt. 2 (1990) 253-265 



2 



1. Introduction 

The generalised coherent states introduced by Titulaer and Glauber and discussed 
by Stoler [2], differ from coherent states by the extra phase factors appearing in the 
decomposition of such states into a superposition of Fock states. Bialynicka-Birula |3] 
has shown that, under appropriate periodic conditions, generahsed coherent states 
go over into discrete superpositions of coherent states. She also has shown how to 
calculate the coefficients of such a superposition. Recently, Yurke and Stoler [1], 
and Tombesi and Mecozzi |5j have discussed the possibility of generating quantum 
mechanical superpositions of macroscopically-distinguishable states in the course of the 
evolution of the anharmonic oscillator. The anharmonic oscillator model was earlier 
used by Tanas [6] to show a high degree of squeezing for large numbers of photons. The 
two-mode version of the model was used by Tanas and Kielich [7j to describe nonlinear 
propagation of light in a Kerr medium, predicting a high degree of what was called 'self- 
squeezing' of strong light. The comparison of quantum and classical Liouville dynamics 
of the anharmonic oscillator was made by Milburn [8], and Milburn and Holmes [9]. 
Kitagawa and Yamamoto [lOj have used the model in their discussion of the number 
phase minimum uncertainty state that can be obtained in a nonlinear Mach-Zehnder 
interferometer with a Kerr medium. They introduced the name 'crescent squeezing' 
for the squeezing obtained in the model, to distinguish it from 'elliptic squeezing' of an 
'ordinary' squeezed state. The terms 'crescent' and 'elliptic' stem from the shapes of the 
corresponding contours of the quasiprobability distribution Q{a, a*, t). The anharmonic 
oscillator model has also been discussed by Pefinova and Luks [llj from the point of 
view of photon statistics and squeezing. Quantum field superpositions have recently 
been discussed by Kennedy and Drummond [I2], and by Sanders [13]. Some properties 
of generalised coherent states have been discussed by Vourdas and Bishop [Hj. 

In this paper explicit analytical expressions describing superpositions of up to four 
coherent states are obtained for the anharmonic oscillator model with two different 
orderings of operators in the interaction Hamiltonian. The maximum number of clearly 
distinguishable coherent states in a superposition is estimated, and the rules describing 
the sequence of particular superposition states as time elapses are given. The results are 
illustrated graphically for the coherent initial states with the mean number of photons 
equal to 4 or 16, for which the evolution of the quasiprobability distribution (QPD) 
Q{a,a*,t) is used to visualise the formation of superposition states. 

2. The anharmonic oscillator model and its evolution 



The anharmonic oscillator model that we discuss in this paper, is defined by the 
Hamiltonian 

H = huja^a + Hc {C = N, S) (1) 
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where a (a^) is the annihilation (creation) operator and He describes the chosen versions 
of the nonlinear interaction Hamiltonian, which are: 

Hn = -hK^a^Ya^ = -hKh{h — 1) (2a) 

1 1 

Hs = -hK{a^ay = -hKh'^. {2b) 

Here n = a^d is the number of photons operator and k is the coupling constant, which 
is real and can be related to the nonlinear susceptibility x^^^ the medium if the 
anharmonic oscillator is used to describe the propagation of laser light in a nonlinear 
Kerr medium. Both versions of the interaction Hamiltonian are in use, depending on 
the authors. The difference between them seems to be trivial because it means a change 
in the free oscillator frequency of the oscillator. When the homodyne detection of 
squeezing is to be applied, however, this extra phase shift can be significant in the long- 
time limit [15]. Thus, the question arises: which version is to be used in a particular 
physical situation? From the point of view of quantum-classical correspondence the 
normal ordering is preferable because it makes the transition from the quantum to the 
classical description via coherent states quite transparent and preserves the classical 
meaning of the nonlinear susceptibility of the medium. However, since both versions 
are used in the literature, in this paper we consider both of them separately in order to 
make the difference more explicit. 

Since the number of photons n is a constant of motion (it commutes with both 
versions of the interaction Hamiltonian) the state evolution of the system is described, 
in the interaction picture, by the Schrodinger equation 

ihfuc{t) = HcUc{t) (3) 

where Uc{t) is the time evolution operator. In the propagation problem of light 
propagating in a Kerr medium, one can make the replacement t = —z/v to describe the 
spatial evolution of the field, instead of the time evolution. The solution of equation ([3]) 
is given by 

lJc{z) = exp i^Hc) = exp (idc{n)) (4) 



where 



or 



and 



T 

^N{n) = -h{n 



9s{h) = Jn^ (5) 



2 



T = KZ/V (6) 

is a dimensionless length of the medium (or time in the time domain). Since the 
difference between the time and the spatial descriptions is trivial (the main effect is 
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the change in sign), we write formulae for the spatial description but use the terms time 
or length interchangeably in the text. 

If the state of the incoming beam is a coherent state lao), the resulting state of the 
outgoing beam is given by 

mr)) = Uc{r)\ao). (7) 
Using the well known decomposition of the coherent state |ao), 



oo 



«o)=exp -^ E^H (8) 



2 / n=o VriJ. 
we obtain from equations and 



oo 



where 



(r)) = |ao, r) = exp ( ) I] ^ ^xp [iedn)] \n) (9) 

eN{n) = ''-n{n-\), Osin) = (10) 

Because of the presence of the additional phases O^iji) or 9s{n), the resulting state is 
a generalised coherent state [H [2] which can be, under certain conditions [3] , a discrete 
superposition of coherent states. Some of these superpositions will be given in the next 
section. 

A good representation of the field state resulting during the evolution of the 
anharmonic oscillator is the quasiprobability distribution Q{a,a*,t) defined as, see 0: 

Q{a,a*,t) = Tr(p(r)|a)(a|) = (a|p(r)|a). (11) 

This function satisfies the relations 

d^a 

Q{a,a*,t) = 1 (12) 

TT 

and 

< g(a*,a,t) < 1. (13) 

The properties of this function for the anharmonic oscillator both in the classical and the 
quantum descriptions of the oscillator have been discussed by Milburn [S] and Milburn 
and Holmes 

In the case of the initial state |ao), the Q-function has the form 

Q(a, a*, 0) = exp(— |a — aoP) (14) 

which is a Gaussian bell centred on oq. 

Since the state of the outgoing field is given by equation ([9]), its density operator is 
p = \iP{t)){iP{t)\, and the corresponding quasiprobability distribution is given by [SjlTU]: 

Qcia,a*,T) = {a\^lj{T)){'ip{T)\a) 



n=0 



exp(-|ap - laoHl ^ ^^P ['^^c{n)] (15) 



2 
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where 6c are given by equations (fTU]) . This quasiprobabihty distribution will be 
illustrated graphically for some specific values of r to show the formation of the 
superposition states. 

According to equations ([9]) and (fT5|) . it is clear that both the state itself and the 
quasiprobabihty distribution exhibit periodic behaviour, however, the two versions of 
the anharmonic oscillator have different periods. Since n{n—l) is always an even number 
(contrary to n^, which can be odd) we have from ffTOj) that the period for the normally 
ordered version of the interaction Hamiltonian is one half of the period for the 'squared' 
version. We have 

Qc{oi,(y* + T) = Qc{a,a* ,t) (16) 
with the periods 

Tjv = 27r, Ts = 47r (17) 

The periodic behaviour of the quasiprobabihty distribution ( fTSl) . or the state ([9]), 
can be observed in the long-time (or long- length) limit. Estimates based on realistic 
values of the non-linear susceptibility of the medium give, for a length of the medium 
of the order of one metre, values of r of the order of 10~^ |7j. This makes it rather 
unrealistic to observe periodic behaviour, at least in the case of the Kerr medium. Such 
a periodic behaviour is, on the other hand, an essential feature of the quantum dynamics 
of the system and, thus, worth studying in its own right. Some quantum features of the 
system such as squeezing are more likely to be observed for a large number of photons 
in the short-time limit [6l [7]. The generation of superpositions of the macroscopically- 
distinguishable states, however, needs rather long evolution times. 

3. Generation of discrete superpositions of coherent states 

The state of the field obtained as a result of the evolution of the anharmonic oscillator, 
which is given by equation ([9]), is a generalised coherent state [H |2]. When the phases 
Ociji) satisfy periodic conditions 

exp [iOcin + A^)] = exp [iOc{n)\ (18) 

for every n, with N being an arbitrary positive integer number, the state ([9]) can be 
represented as a discrete superposition of coherent states |3] 

N 

iV'(r)) = |ao, t) = Y^ Cfcl exp(iv5fc)ao) (19) 

k=l 

where the phases ipk and the coefficients Ck are to be found. For a specific choice of the 
time (length) 

Ts = Ts/N = An/N (20) 

the periodic conditions ( fTSj) are satisfied and, for A^ odd, the superposition ( fT9l) can be 
found directly according to the formulae given by 

ipk = 27ik/N k = l,2,...,N Nodd, (21) 
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and the coefficients Ck can be found from the following system of equations: 



N 



Ckexp{inipk) = exp [i6s{n)] = exp i{2'n/N) 



n 



n 



0, 1,...,A^- 1. (22) 



k=l 



With the choice f l20l) of the time (length), there is some additional symmetry that 
allows us to reduce the number of equations that are needed for finding the coefficients 
Cfe. Using the relations 

exp [ies{N - n)] = exp [iOsin)] (23) 

one easily finds that 

CN-k = Ck (24) 

which means that the number of equations is reduced to ^{N — 1) + 1. In the case of odd 
N, we obtain a superposition of coherent states with their a satisfying the relation 



"0 



When is even, the following relations hold: 

exp [iOsin + N/2)] = (-1)^/^ [i0g(n)] 

exp [ies{N/2 - n)] = (-1)^/^ exp [tOsin)] . 



(25) 
(26) 



If the relation (!25|) is applied to equation (l22l) . it becomes evident that, depending 
on whether A^/2 is odd or even, only the coefficients with odd or even k survive. This 
reduces the number of equations by one half. Instead of equations that are needed 
in the general case, there are only A^/2 equations and, correspondingly, the resulting 
superposition has only A^/2 coherent states. In effect we obtain: 

(i) For A^ even and A^/2 odd 

(^fc = 27r(2A; - 1)/A^ A; = 1, 2, A^/2 (27) 



N/2 

C2fc_i exp [my^fc] = exp [i6s{n)] = exp i{27i/N)n^ 

k=l 



n = 0,l,...,^Ar-l(28) 



and, when the relation ( ]26i) is exploited, the following relation between the coefficients 
is found: 



CAr-(2fc-l) — C2k-1 

which reduces the number of equations to |(|A^ — 1) + 1. 
(ii) For A^ even and A^/2 even 

^k = 27i{2k/N) k = 1,2,..., N/2 

N/2 

C2k exp [impk] = exp [iOs{n)] = exp i{2'K/N)ii? 

k=l 



and again the relation fl26|l leads to 

CN-2k = C2k 



(29) 

(30) 
(31) 

(32) 
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reducing the number of equations to + 1. 

Of course, the numbering of the coefficients C2k-i and C2k can be replaced by c^. 
However, we keep the above notation for Ck in order to indicate their origin and to show 
clearly their symmetry. 

It is evident from the results obtained above that the symmetry of the system under 
consideration plays a crucial role in reducing the problem of finding the superposition 
states. It is also clear from fl2T|) to f l27|) that superpositions of say three states appear 
both for = 3 and for = 6. These are, however, different states. Comparison of ( |2T|) 
and ( 127|) shows that the phases of the two superpositions differ by tt/3, which means 
refiection with respect to the Ima axis. In fact, if we take the time equal to 2Ts/6 = Ts/3 
in f l28|) . we easily recover the state obtained from f l22|) for Ts/3, if we simultaneously 
replace ipk by 2ipk- So, the number of components in the superposition depends on what 
fraction of the period we take for the time. If the fraction of the period is m/N, assuming 
that this fraction cannot be reduced, the number of components is equal to N for N odd, 
and N/2 for N even. If the fraction m/N can be reduced, the above rules must be applied 
to the reduced fraction. It is also not difficult to prove that coefficients Cfc obtained for 
(A^ — k)/N are complex conjugates of those for k/N. These are general rules governing 
the process of generation of the discrete superpositions of coherent states during the 
evolution of the anharmonic oscillator. Of course, when the evolution starts at time 
r = 0, the superpositions with large numbers of components will appear first, while the 
superpositions with few components cannot appear before the time approaches half the 
period or a fraction of this with small denominators (2,3,4,...). Since the superpositions 
with a small number of coherent components are most interesting in the discussion 
of macroscopically-distinguishable states, we give here some examples of such states 
obtained with the use of the formulae derived in this section. 

For N = 2, according to the rules, there is only one coherent state in the 
superposition, which according to f l271) and ( l28l) is equal to 

|ao, Ts/2)s = \ao, 2n)s = \ exp{in)ao) = | - ao). (33) 

When = 4, the number of components is equal to two and, according to ( l30i) and ( l3Ti) 
we have 

\ao,Ts/4:)s = \ao,n)s = [exp(-z7r/4)| - ao) + exp(i7r/4)|ao)] • (34) 

Equations (133|) and ( |34|) are the results given by Yurke and Stoler |1] in their discussion 
of the problem of generation of a superposition of macroscopically-distinguishable states. 
For r = 3T5/4 we obtain the state with coefficients that are complex conjugates of the 
coefficients in flM|) . 

For A^ = 3, we have obtained from (^I^ and fl2^ the following superposition: 
100,^5/3)5= |ao,47r/3)5= --^[exp(-i7r/6)| exp(i27r/3)|ao) 

+ exp{—in / 6) \ exp(— ■j27r/3)ao + ^l^o)] (35) 
and for r = 2Ts/3 the state with the complex-conjugated coefficients is obtained. 
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For = 6, according to (1271) and (I28p . we have 
\aQ,Ts/6)s = \ao,2n/3)s = --^[exp(z7r/6)| exp(ii7r/3)|ao) 

+ exp(i7r/6)| exp(— i7r/3)ao) ~ "^1 ~ <^o)]- 



(36) 



This state is different from (135!) . but the state (135|) is obtained for r = 2T5/6 = T5/3. 
For = 8, we have obtained: 



ao, 3X5/8)5 = |ao,37r/2)5 

= - Wiao) + exp(-z7r/4)| - ao) + I - ^"o) - exp(-i7r/4)|ao)] • (38) 



The results (137|) and (!38|) are interesting because they correspond to the plots of the 
QPD Q{a,a*,t) given, for = 2.0, by Milburn [8], from which four Gaussian peaks 
of the QPD are clearly visible. The two-peak structure corresponding to the state (|35|) 
is also evident. Knowing the superposition states makes the interpretation of the 
multipeak structure of the QPD quite transparent. The quasiprobability distribution 
has four peaks because it represents a superposition state composed of four coherent 
states. This, of course, is true when the component states are well separated and the 
interference terms are negligible. Thus, the question arises: when can the components of 
the superposition be considered as well separated? To answer this question we have to 
remember that the Gaussian quasiprobability distribution (IT^ representing a coherent 
state has a finite width. If we assume, somewhat arbitrarily, that the states are well 
separated when the distance between their Gaussian peaks in the complex a-plane is 
equal to the diameter of the contour obtained when the section of the Gaussian bell is 
made at 0.1 of its height, the diameter is then estimated by the value 2(ln 10) = 3.03. 
On the other hand, all the coherent states entering the superposition have their 
parameters such that \ak\ = |ao|- This means that all the Gaussian peaks representing 
such states are distributed regularly around a circle of radius |q;o| in the complex a- 
plane. So, the maximum number A^max of the well-separated Gaussians for given |q;o| 
can be estimated by 



This estimation gives for |q;o| =2 that the maximum number of well-separated peaks 
in the QPD Q{a,a*,t) is four. These are the four peaks obtained by Milburn [8]. In 
fact, the five-peak or even six-peak structure of the QPD can still be identified when 
the proper r is taken, but the peaks are not well separated and their shapes are strongly 
affected by the interference terms. 

To make this point clearer we write down the analytical expression for the QPD 
of the discrete superposition of coherent states which can be split into the sum of pure 
Gaussians and another sum describing the interference terms 



(yo,Ts/8)s 



ao,7r/2)5 

- [\iao) - exp(OT/4)| - ao) + I - ^"o) + exp(OT/4)|ao)] (37) 




(39) 




(40) 
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Re a 




Figure 1. Plots of the quasiprobability distributed Q{a,a* ,t); (a) for ao = 2.0 and 
T = Tn/A (the four-peak structure) and (b) for ao = 4.0 and r — Tjv/8 (the eight-peak 
structure). 



where 



with 



N 



Qcanss ^ ^ Qk 
k=l 

Qint = X! ^ Qkl 

k>l 



Qk = Icfcpexp - 



(41) 



(42) 



2ReQki = 2|cfc||Q| exp [ --\a - ak\ 

X cos [7fc - 7/ + |a||ao|(sin A^jfc - sin A(^/)]) 



-\a — aA 



In equation f H5]) we have used the notation 

Cfc = |cfe| exp(z7fc), OL = \a\ exp{i(p), ao = \ao\ exp{iipo) 

and 



(43) 



(44) 



Aifik = fk - fo - f- (45) 

In deriving equations the superposition state (HM has been used. It is 

clear from fH3|l that due to oscillations of the cosine function the interference terms can 
have a number of peaks. However, due to the exponential factor the amplitudes of these 
peaks are very small whenever the states k and / are well separated. 'Well separated' 
means here that la^ — a/p ^ 1. This gives us another criterion for good separation 
of states. In the following we illustrate the formation of the superposition states by 
showing pictures of their QPD for special situations. 

Before doing this, however, we have to make some comments on the behaviour of the 
normally ordered version of the anharmonic oscillator. Our analytical formulae ( 133|1 - 
( l38l) are for the 'squared' version of the anharmonic oscillator; however, it is clear 
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Figure 2. Contours of the quasiprobability distribution for the 'squared' version of the 
anharmonic oscillator obtained for sections at j, ^ and | of the height. The parameters 
are ao = 2.0, and (a) r equal to Ts/8, (b) Ts/6, (c) Ts/4, (d) Ts/S, (e) 3Ts/8 and (f) 
Ts/2. 



from (fTOj) that ONiji) = Osiji) — irn. If this is inserted into equation (j9]), it is seen 
that the superposition states obtained from the normally ordered version acquire an 
additional phase ip = — r/2. Thus, the only change that is needed to obtain the results 
for this version is the replacement of the ipk by the — |r in all the formulae obtained 
in this section. Geometrically this means the rotation of the QPD picture by the angle 
if = —t/2 in the complex a-plane, without any change in its shape. This rotation may, 
nevertheless, change in an essential way the general view of the QPD. For example, the 
QPD representing the state fl33l) is a Gaussian centred at — ao, and after the rotation 
by vr it becomes a Gaussian centred at +ao, that is at the initial value of ao. There 
is no single Gaussian peak at —ao for the normally ordered version. This is a general 
rule, related to the fact that the period for the normally ordered version is one half of 
the period for the 'squared' version; this will be convincingly shown in the figures. 

In figure 1 we show two examples of the QPD. The first example exhibits four 
Gaussian peaks for the case considered by Milburn of ao = 2.0. The second example 
shows the eight Gaussian peaks in the case ao = 4.0, which according to our 
estimate (139|1 is the maximum number of well-separated Gaussians in this case. Both 
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Figure 3. The same as in figure 2 but for the normally ordered version of the 
anharmonic oscillator. To visualise the differences between the two versions of the 
model, the same moments, r, arc taken, that is (a) r = T5/8, (b) Ts/6, (c) Ts/'i, (d) 
Ts/3, (e) 3Ts/8 and (f) Ts/2. 




Figure 4. The same sections of the QPD for the normally ordered version of the 
anharmonic oscillator as in figure 3, but for the parameters: ao = 4.0 and (a) r equal 
to Tn/8, (b) Tn/6, (c) Tn/5, and (d) Tjv/S. 
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examples represent quite regular shapes confirming the formation of superposition states 
with a definite number of components. All examples presented in our pictures are 
obtained numerically from the expression ( |T5i) . In figure 2 contours of the sections at 
|, |, I of the height of the QPD for the 'squared' version of the anharmonic oscillator are 
presented, for ao = 2.0, showing that for the maximum number of well-separated states 
Nuiax = 4 [estimated according to fl39l) ] the contours are not very regular circles yet, as 
they should be for the independent Gaussians. However, the regular four-peak structure 
is clearly visible. As time elapses the structures with various numbers of peaks appear 
which represent the superposition states given by the formulae (l33l) - (!38l) . The smaller 
the number of peaks the better is the separation of the states, and the more regular is the 
QPD. The sequence of the pictures is obtained for r = T5/8, Ts/6, Ts/4:, Ts/3, 3Ts/8 
and T5/2, respectively. The identification of the QPD structures with the corresponding 
superposition states is quite obvious. In figure 3 the same sequence of the QPD 
structures is shown for the normally ordered version of the anharmonic oscillator. The 
differences between the two versions are quite evident. It is seen that the period for the 
normally ordered version is really one half of the period for the 'squared' version. We 
have chosen the same sequence of r for both versions to visualise the differences. The 
structures obtained for the 'squared' version after half the period, which is T5/2 = 2tt, 
become rotations of the initial structures by the angle of vr, while for the normally 
ordered version after the same time the initial structures are recovered. In figure 4 
the contours of the QPD are presented for ao = 4.0 and the normally ordered version. 
Again the structure with the maximum number of well-separated states, which in this 
case is iVmax = 8, shows some irregularity, but as the number of peaks decreases 
the structures become more and more regular. The pictures have been obtained for 
T = Tjv/8, Tjv/6, T/v/5 and Tjv/3, respectively. We have chosen here, as examples, two 
structures with even number of peaks and two structures with odd number of peaks, 
only. Figures 2-4 are on the same scale, which shows that the radius of the circle around 
which peaks are located is in figure 4 twice as large as that in figures 2 and 3, whereas 
the radii of the individual Gaussian bells are the same. One should also remember that 
the coefficients of a superposition with N peaks scale as \ck\ = 1/\/N. This means 
that the heights of the peaks are N times lower from the single coherent-state Gaussian. 
If the initial number of photons |aoP becomes large, |ao| is large, and the maximum 
number of well-separated states iVmax is also large, but the amplitudes Ck of these states 
become smaller. 

4. Conclusions 

In this paper we have analysed the process of the generation of discrete superpositions 
of coherent states in the course of the evolution of the anharmonic oscillator. Two 
versions of an anharmonic oscillator that are used in the literature have been compared 
from the point of view of forming the superposition states. It has been shown that under 
an appropriate choice of the evolution time as a fraction of the period, the symmetry 
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inherent in the system permits a considerable simphfication of the problem of finding 
the coefficients of the superpositions. The number of equations that must be solved 
is drastically reduced by the symmetry. Some examples of the superposition states 
have been obtained analytically for superpositions of up to four components. General 
rules governing the formation of superpositions with a definite number of states have 
been given. The process of the formation of superposition states has been illustrated 
graphically by showing pictures of the corresponding quasiprobability distributions. The 
regular structure of the QPD that is obtained when the time (or length) becomes a 
fraction of the period can be easily interpreted as representing the superposition state 
that occurs for this time. The maximum number of well-separated states for given 
|q;o| has been estimated. This number becomes large when |ao| becomes large, and 
regular structures of the quasiprobability distribution with a large number of Gaussian 
peaks can be obtained. Some of these structures have been shown in the figures. 
Structures with both even and odd numbers of peaks are possible. Our results shed some 
new light on the problem of generating discrete superpositions of coherent states and 
make a contribution to the discussion of the possibility of generating macroscopically- 
distinguishable quantum states [H E] as well as to the problem of the phase space 
interferences [HI [17] . 
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